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Abstract 

Theoretical predictions for B decay rates are rewritten in terms of the T(IS') 
meson mass instead of the b quark mass, using a modified perturbation ex- 
pansion. The theoretical consistency of this expansion is shown both at low 
and high orders. Our method improves the behavior of the perturbation se- 
ries for semileptonic and nonleptonic inclusive decay modes, as well as for 
exclusive decay form factors. The results are applied to the determination of 
the semileptonic B branching ratio, charm counting, the ratio of B — * Xtv 
and B — > Xev decay rates, and form factor ratios in B — * D*ev decay. We 
also comment on why it is not possible to separate perturbative and nonper- 
turbative effects in QCD. 
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I. INTRODUCTION 



Heavy hadron decays can be computed using a systematic expansion in powers of a s (mq) 
and AQ CD /mQ, where mq is the mass of the heavy quark and Aqcd is the nonperturbative 
scale parameter of the strong interactions. In the mq — ► oo limit, inclusive decay rates are 
given by free quark decay and the order Aqcd /^iq corrections vanish The leading 

nonperturbative corrections of order AQ CD /mQ are parameterized by two hadronic matrix 
elements 0-§[ Ax and A2. These results are now used to determine the CKM matrix elements 
\Vcb\ and \Vub\, using experimental data on inclusive semileptonic B meson decays. For 
exclusive decays at lowest order in I/tuq, one can compute form factors in terms of a 
perturbative series in a s times the Isgur-Wise function [0 . Nonperturbative corrections are 
parameterized by hadronic matrix elements, and are of order (Aqcd/^q)™- In some cases 
at zero recoil, the order l/mq correction vanishes ||. 

At present, the largest theoretical uncertainties in the computation of inclusive B — > 
X c eu and B — ► X u ev decay rates arise from poor knowledge of quark masses. The pole 
mass of a heavy quark is an infrared (IR) sensitive quantity, and is not well defined beyond 
perturbation theory [§,0. This is related to the bad behavior of perturbative corrections 
to the inclusive decay rate when it is written in terms of the pole mass |p| , |12| . A simple 
method of avoiding problems with the quark mass is to use instead the B mesons mass. Un- 
fortunately, the B meson and b quark masses differ by order Aqcd, and so this reintroduces 
a Aqcd/ 7 ™?) correction to the inclusive decay rate. The decay rate has also been rewritten in 
terms of the infrared safe MS mass of the b quark, but the uncertainties still remain sizable 
because the value of the MS mass is not known very accurately, and because the first few 
terms in the perturbation series still remain large. 

In this paper, the theoretical predictions for semileptonic B decays are rewritten in terms 
of the T(15*) meson mass rather than the b quark mass. At the same time a modified per- 
turbation expansion (referred to as the upsilon expansion) must be used, which is discussed 
in detail in Sec. II. This procedure eliminates the uncertainty due to the mf factor in the 
decay rates, and at the same time improves the behavior of the perturbation series. Our 
formulae relate measurable quantities to one another and the resulting perturbation series 
is free of renormalon ambiguities. We will also see numerically that the perturbative correc- 
tions are small when the B decay rate is written in terms of the T mass. We have discussed 
the procedure in an earlier publication |13 |, where it was applied to inclusive semileptonic 



B decays. In this article, we describe the method in more detail and also explain why it 
leads to an improvement in the perturbation series expansion. The method is applied to 
inclusive nonleptonic decays as well as to exclusive decays. The perturbation series is better 
behaved for all the processes we have examined. The semileptonic and nonleptonic decay 
computations can be combined to obtain new values for the semileptonic branching ratio 
and the average charm multiplicity in B meson decay. 

The outline of this paper is as follows. In Section II we define the upsilon expansion. In 
Sec. Ill we show that it is the theoretically consistent way of eliminating the b quark mass 
from the combination of the perturbation theory results for B decay rates and the T mass. 
We also comment on the role of nonperturbative effects in the theoretical expression for the 
T mass. Exclusive decays are discussed in IV. In Sec. V we give the results for inclusive 
semileptonic and nonleptonic decay rates in the upsilon expansion, the ratio of decay rates 
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R T = BiB — > X c rv)jB[B — > X c ez/), the semileptonic B branching ratio, and charm counting 
in B decay. Sec. VI contains our conclusions and comments on further applications of our 
results. 



II. THE UPSILON EXPANSION]] 

The results for the OPE calculation of inclusive B decay rates can be schematically 
written as 



ijk) 



Gp |Vckm| 

1927T 3 



ml {PS) 



1 + — e + 



+ 



7T 
1 

mi 



7T 



(A) 



(Ai + A 2 ) + . . . 



7T 



e* + ... 



(1) 



where the precise coefficients are not shown, ijk denotes an arbitrary final state created by 
the b — > ijk transition at short distance (e.g., ijk = civ, ccs, etc.), Vckm contains CKM 
matrix elements, PS is the b quark decay phase space at tree level including color factors 
(e.g., ijk = uud is a factor 3 enhanced compared to ijk = ueu), m b is the b quark pole mass, 
f3 = 11 — 2nf/3 is the first coefficient of the QCD /^-function, and a s is the running coupling 
constant in the MS scheme at the scale /i. It is convenient to expand the coefficient of a 
given a s order in powers of f3$ which is related to the Brodsky-Lepage-Mackenzie (BLM) 
prescription For nonleptonic and rare B decays there is another large parameter in 

Eq. (|1|), \n(m w / m b ) ~ 2.8, which will be discussed in Sec. IV. These large logarithms can be 
summed using the renormalization group equations. The reason for introducing the variable 
e = 1 is explained in the next paragraph. 

Schematically, the perturbative expansion of the T mass in terms of m^ is 
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where £ = ln[/x/ (mbCt s Cp)], Cp = 4/3, and the precise coefficients are again not shown. Note 
that this series is of the form {a 2 s1 a^/3 , a^Q, ■ ■ .}, whereas the corrections in Eq. (|1|) are 
of order {a s , a 2 p , cx^Pq, ■ ■ ■}■ We will show in Sec. Ill that to combine these two equations 



B decays written in terms of the T mass should be referred to as upsilon-expanded, even though 
the expansion parameter will be denoted by e. Blame this on Donald Knuth (or on the Greeks?), 
as v ($\upsilon$) is similar to v ($v$) in TgX. 

2 The subscript BLM will be used to denote terms of order e n /?Q _1 in the perturbation expansion. 
In some cases, the entire e 2 contribution has not been computed, and only the BLM piece is known. 
In cases where the entire e 2 term is known, the BLM contribution is about ten times larger than 
the non-BLM part. 
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in a theoretically consistent manner, terms of order a" in Eq. ([J) should be viewed as if 
they were only of order a™ -1 . For this reason, the power of e in Eq. @ is one less than 
the power of a s . It is also convenient to choose the same renormalization scale, //, when 
Eqs. (H) and (@) are combined. The prescription of counting [a s (mj]' 1 in the B decay rate 
as order e n , and [a s (m(,)] n in the T mass as order e n_1 will be called the upsilon expansion. 
It is expected that the infrared sensitivity present separately in Eqs. (|l|) and @ will cancel 
to all orders in perturbation theory in e. Some arguments that support this will be given in 
the next section. In the upsilon expansion Eq. (p]) takes the form 
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where PS' means that the b decay phase space is also expanded in e. We emphasize that 
the coefficients in the series in e contain different orders in a s . We show in Sec. Ill that this 
is theoretically consistent. Our results in Sec. IV and V demonstrate that this expansion 
simultaneously eliminates the uncertainty in the theoretical predictions related to the depen- 
dence on the poorly known quark masses and also improves the behavior of the perturbation 
series. 

The precise form of Eq. (|lj) depends on the decay channel under consideration. The 
expression for the T mass in terms of mj, is given by |15|,|l6| , 
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The ellipsis denote terms of order a^, and non-BLM terms of order a A s . For /i of order mj, 
Eq. (^) shows no sign of convergence; for /i = m6 it yields mx = 2m?, (1 — 0.0 lie — 0.0 16e 2 — 
0.024blm£ 3 — • • •)> using a s (m;,) = 0.22. (This value of a^mj) will be used whenever we 
present numerical results.) The bad behavior of this series is unimportant (it would be 
somewhat better in terms of the MS mass at a low scale). The only physical question is 
what happens when we use Eq. (H) to predict B decay rates in terms of mx- 

Finally, decays with a charm quark in the final state depend both on and m c . It is 
convenient to express these decay rates in terms of mx and Ai instead of m& and m c , using 
Eq. (D and 

m b -m c = m B -m D + f-^ 7^=^) + ■■■ , (5) 

\2m B 2m D J 

where TUb = (3ms* +ms)/4 = 5.313 GeV and % = (3m_D* +"1d)/4 = 1.973 GeV. One 
could imagine replacing m c by the J/ip mass using a procedure similar to that for the b 
quark. However, the nonperturbative errors in this case are not under control. 



III. THEORETICAL CONSISTENCY 



In Ref. |T3| we showed that our method of combining the expansion of inclusive B decay 
rates with the expansion of the T mass to eliminate the dependence on the b quark mass 
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is consistent at large orders in perturbation theory in the BLM approximation. In this 
section we briefly repeat this argument (Sec. A) and show that the upsilon expansion is 
also consistent at first and second order in e, including non-Abelian contributions (Sec. B). 
Our arguments prove that, as far as the cancellation of infrared sensitive contributions in 
the perturbative expansion of the B decay is concerned, the upsilon expansion is equivalent 
to the use of a short-distance mass like the MS mass. They also explain why the upsilon 
expansion requires powers of a s to be counted differently in mx and B decay rates. We also 
comment on the relation between perturbative and nonperturbative effects in the T mass 
(Sec. C). 



A. Renormalons in the large [3q approximation 



The pole mass m b is an infrared sensitive quantity, which is not well defined beyond 
perturbation theory. It can be related to an infrared safe mass, such as the MS mass via 
(for n f = 4) fl7[ 
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This relation has terms of the form o.™/3q 1 n\ at high orders, leading to a renormalon ambi- 
guity. The Borel transform of the perturbation series relating the pole mass to the MS mass 
is 1,0 
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where C is a scheme dependent constant (C = 
ambiguity is given by the u = 1/2 pole, yielding 



-5/3 in MS). The leading renormalon 
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The a s perturbation series for the inclusive decay rate written in terms of the pole mass 
also has a renormalon singularity. As a result, the inclusive decay perturbation series is 
badly behaved at high orders. Renormalons cancel when the inclusive decay rate is written 
in terms of the MS mass ||l8|-f2l|1 , rather than the pole mass. While this cancellation is present 
at high orders, the perturbation series written in terms of the MS mass at a renormalization 
scale fj, of order m b still contains large corrections at low orders (see, Ref. |TT| , p!2| and some 



examples in Sec. V). On the other hand, the perturbation series written in terms of the MS 
mass at a low scale has a better behavior. But such a choice still does not remove the quark 
mass uncertainty in the decay rates. 

The expansion of the T mass also becomes better behaved when written in terms of the 
MS mass at a low scale (of order the Bohr radius, ma s ), rather than in terms of the pole 
mass. This fact is important for extracting the short distance masses of the bottom and top 
quarks from the T system and from the location of the peak of the IS ti resonance in the total 
cross section for ti production in e + e~ annihilation close to threshold. One might think that 
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the numerical results obtained in the upsilon expansion are equivalent to the result obtained 
by extracting a short-distance b quark mass from the T mass and then using it to predict 
B decay rates. Due to the fact that only the first few terms in the perturbative expansions 
are known, this is not the case because in this procedure the information on the correlations 
between the short-distance mass, the strong coupling constant and the renormalization scale 
is lost. The upsilon expansion, which completely eliminates the intermediate step of using 
any quark mass takes full advantage of these correlations and leads to smaller uncertainties. 

The potential between two static quarks in a color singlet state, V{r) = —Cpg^/r, also 
has renormalon ambiguities. The Borel transform is given by 



The leading renormalon ambiguity of V also comes from the u = 1/2 pole, 

Al/ = ^e- c % CD . (10) 
Po 



It follows from Eqs. (|7|) and (0) that the order Aqcd ambiguity cancels in 2m?, + V p3|j2^ 

Equation (|9]) also explains the mismatch between the power of a s and e in Eq. (|j), which 
may seem somewhat unnatural at first sight. The coefficient of a" +2 /3g in Eq. ([2]) is given 
by the T matrix element of the coefficient of u n in the expansion of V(u) in Eq. (|9|) about 
u = 0, multiplied by n\. For large enough n, up to corrections suppressed by at least one 
power of n, the terms in this expansion of V(r, u) are proportional to 

(2tiT Sl^f ~ p (2„)». (11) 
r p=0 P- 

This explains why in the BLM approximation and in higher orders in e, the terms in Eq. (0) 
are generically of the formf] 

/ pn pn-1 \ 

(12) 

since the T matrix element of ln(/ir) ~ £. For large n, 

+ ... + l)^exp(£) = —V—, (13) 
J m b a s C F 

reducing the difference between the powers of a s and /?o to one, thus eliminating the mis- 
match in the large order behaviors of Eqs. ([I]) and (fj). This has to happen since my is a 




3 The actual coefficients of the logarithms in Eq. (g) differ from those displayed in Eq. (12) 
because of multiple insertions of the corrections to the pure Coulomb potential, —Cpa s jr. The 
contributions arising from multiple insertions, however, only constitute corrections suppressed by 
at least one power of n. 
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physical quantity, so the renormalon ambiguities must cancel in Eq. (|2]) between 2rrib and 



the potential plus kinetic energies p3jp4 



If NRQCD exists as a well defined low energy effective field theory, then the AQ CD r 2 
renormalon in the potential corresponding to the u = 3/2 pole in Eq. (H) must cancel 
against the leading renormalon in the kinetic energy contribution to the T mass. While this 
has not been shown explicitly, it is quite conceivable that this happens. Studying higher 
order renormalons is less interesting as there are certainly nonperturbative contributions of 
order AQ CD /(ma s ) 3 to m T (e.g., from the vacuum expectation value of G 2 ^ p5|-p7|), and so 
higher order renormalon effects are of the same order as other nonperturbative effects which 
are known to exist. 



B. Infrared sensitivity at low orders in e 

One can study the infrared sensitivity of Feynman diagrams by introducing a fictitious 
infrared cutoff A of order Aqcd- The infrared sensitive terms are nonanalytic in A 2 , such as 
(A 2 )"/ 2 or A 2 " In A 2 . These terms arise from the low-momentum part of Feynman diagrams. 
This is similar to what occurs in chiral perturbation theory, where the nonanalytic terms 
in the quark mass (or equivalently, the pion mass-squared) come from the low-momentum 
behavior of Feynman diagrams. Diagrams which are more infrared sensitive, i.e., have 
contributions (A 2 )™/ 2 or A 2 ™ In A 2 for small values of n, should have larger nonperturbative 
contributions. In particular, linear infrared sensitivity, i.e., terms of order A, are a signal of 
Aqcd effects, quadratic sensitivity, i.e., terms of order A 2 In A 2 are a signal of Aq CD effects, 
etc. 

In the following we show that the linear infrared sensitivity present separately in the 
expansion of the B decay rate and in the expansion of the T mass in terms of the b quark 
mass cancel at least to order e 2 when they are combined using the upsilon expansion. (At 
order e this is equivalent to the cancellation of the leading renormalon ambiguity in the BLM 
approximation, discussed previously.) The argument is essentially a combination of those in 
recent papers by Beneke ||24j| , and by Sinkovics, Akhoury, and Zakharov [21j. As far as the 



cancellation of infrared sensitive contributions in the perturbative expansion of B decay is 
concerned, the upsilon expansion is therefore equivalent to a scheme where a short- distance 
mass like the MS mass is used. 

The proof proceeds in two steps. First one notes that the linear IR sensitivity cancels 
in inclusive decay widths of a heavy quark when it is expressed in terms of an infrared 
safe short distance mass (such as the MS mass). To leading order in a s this was shown in 
Refs. [P^-pD|,|5|] . To order a 2 , where non-Abelian contributions first occur, the cancellation of 



the linear IR sensitivity was shown only recently in Ref . |21| . This cancellation of the linear 
IR sensitivity is expected to hold at higher orders in a s as well, but the demonstration of 
this appears highly non-trivial. The second step is to show the cancellation of the linear IR 
sensitivity when the T mass is expressed in terms of an infrared safe short distance mass. 
This was done by Beneke ||24|| , and we will repeat some of his arguments below. The reason 
is to explain why it is imperative to assign one less power of e to each term in Eq. (|j) than 
the power of a s even at low orders in e. 

Let us first consider the leading contribution to the heavy quark self energy, shown in 
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— £_ — a — » a , 

FIG. 1. One-loop correction to the heavy quark self energy, and the tree- level contribution to 
the QQ potential. 



Fig. Since we are only interested in the infrared behavior, we may calculate the self 
energy diagram in terms of the bare mass, mo, which differs from the MS mass only by 
ultraviolet subtractions. The self energy is given by 

tY _ r n 2 f A + j + ^q)7 m , u s 

F9s J (27iY[{p + q y-ml + ie}{q^ + ie)- 1 1 

To obtain the leading infrared contribution to the mass shift coming from the region of 
small loop momentum, we can set p = m^v, where v is the four-velocity of the heavy quark 
(f 2 = 1). This gives 

5m = —iCp ql [ -, — tt ~, tt^ r ■ (15) 

ys J (2vr) 4 (vq + ie){q 2 + ie) V ; 

The integration over q° is simplest by closing the contour in the upper half plane, yielding 

d 3 q 1 



5m 

2 



(16) 



Of course, this formula is only valid in the infrared. Restricting the range of integration to 
|q| < A, where A is of order Aq C d, it gives 5m ~ a s \. 

Let us next consider the heavy quark potential. It was shown in Ref. that the 
potential in momentum space is free of linear IR sensitivity. The linear IR sensitivity of the 
contribution of the potential to m-f arises in a momentum space calculation from taking its 
matrix element. This occurs because the momentum space wave function ~ 1/(1 + a 2 q 2 ) 2 
does not suppress the low momentum contribution of V(q). In coordinate space, the wave 
function ~ exp(— r/a) is exponentially small at long distances, and the IR sensitivity of the 
contribution of the potential to m-j arises from Fourier transforming V(q) into coordinate 
space to obtain V(r). Therefore, the linear infrared sensitivity of the potential energy can 
be isolated by restricting the range of integration in the Fourier transform 



5V= ( -^U iq - r Wq). (17) 

i|q|<A(27T) 3 y 1 

Choosing A of order Aqcd ; much smaller than 1/a, the e* q r factor can be replaced by 1, and 
the IR sensitive term is a constant times a s A.[] Since this contribution to the potential is 
independent of r, its T matrix element is simply itself, 



J The argument in this and the next paragraph are directly from Beneke [24|. 



5 The assumption Aqcd <C 1/a is essential for the argument to hold. For the T(15) ground state 
this assumption is valid. 
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T 



T)~a s A. (18) 



From Eq. fllq) it is evident that the infrared sensitive contribution in the T matrix element 
of the heavy quark potential is of one lower order in a s than the T matrix element of the 
heavy quark potential itself in the limit A — > 0, 



T 



T)~m 6 af. (19) 



This mismatch arises because the T matrix element of 1/r is of order the inverse Bohr radius 
1/a ~ mf,a s . Consequently, it is the infrared sensitive contribution in the order a 2 s term in 
the perturbative series for my in Eq. (Q) which cancels the infrared sensitive contribution 
of the order a s term in the perturbative series for inclusive B decay rates in Eq. ([!]) (or 
equivalently for the h quark pole mass). It was shown in Ref. |24j] that this cancellation of 
the linear IR sensitivity occurs similarly to all orders in 

Thus it is shown that using the upsilon expansion for B decay rates eliminates all linear 
infrared sensitivity at least to order e 2 . This is the order to which most B decay rates have 
been computed, and the order to which we will present results in Sees. IV and V. While we 
cannot prove that this cancellation is present at all orders in perturbation theory, this is likely 
in view of the fact that it holds both at large orders in the bubble summation approximation 
and in the first two orders including non-Abelian contributions. Consequently, the upsilon 
expansion is insensitive to any order Aqcd ambiguity in the h quark mass which is present 
in perturbation theory. 



C. Nonperturbative corrections to the T mass 

An important theoretical uncertainty in applying the upsilon expansion is the size of 
nonperturbative corrections to the T mass in Eq. (£|). The dynamics of the T system can 
be described using NRQCD ||32|| . The leading nonperturbative corrections to mj arise from 
matrix elements in the T of i/ughtj the Hamiltonian of the light degrees of freedom. In B 
mesons, the leading nonperturbative correction to the B meson mass is due to the matrix 
element of flight) which is the A parameter of order Aqcd- The Aqcd dependence is different 
for the T. f^ght is the integral of a local Hamiltonian density, 

flight = Jd 3 xH hght (x). (20) 

The radius of the T is a ~ l/(m&a; s ), so the matrix element of -flight is of order o 3 Aq CD , by 
dimensional analysis. (Note that the matrix element of -f/iight is of order Aq CD , not m\. Terms 
that grow with m b can be treated using NRQCD perturbation theory.) Using 1/a ~ 1 GeV, 



6 For Beneke's all order argument to hold one has to assume that the heavy quark potential does 
not contain any infrared diverences (or equivalently infrared sensitivity to scales below the inverse 
Bohr radius) to all orders in a s . It was argued in Ref. [28] that this assumption may not hold beyond 
order a^. To order a 3 , i.e. to order e 2 , this assumption is proven by explicit calculations |2(j 
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and Aqcd ~ 350 MeV, of order a constituent quark mass, gives a nonperturbative correction 
of 15 MeV. Using instead Aq CD ~ 500 MeV gives a correction of 60 MeV. We will use 100 MeV 
as a conservative estimate of the nonperturbative contribution to m-f. 

For large enough quark mass, the dominant nonperturbative corrections to the quarkonia 
masses are believed to come from the gluon condensate, (a s G 2 ^ u ). Its contribution to the 



mass of the n 5*-wave state as calculated in Refs. |26|,|27fl is 



, ^ n 6 e n ir(a s G 2 ,) 
(m b C F a s p 

where a s in the denominator should be evaluated at the scale of the order of the Bohr 
radius of the n th state, a s ~ a s (a~ l /n 2 ), and t\ = 1.41, e 2 = 1.59, etc., are numerical 
constants. Values of the gluon condensate in the literature vary in the range (a s G 2 ) ~ 
(0.05 ± 0.03) GeV 4 , and so the uncertainties in Eq. (P) are sizable. In Ref. pf, Eq. (gl]) 
is estimated to give a 60 MeV contribution to the T(IS') mass. Equation (|2T|) seems to 
be enhanced by l/a s compared to our naive estimate of a 3 Aq CD above. This is due to the 
fact that (a s G^) and not (G 2 ) is renormalization group invariant. In our estimate the 
operator 7in g ht is renormalized at the scale fi = a -1 . The renormalization group scaling 
from fi = a" 1 to a low scale fi ~ 1 GeV gives an additional factor a s (lGeV)/a s (a _1 ), which 
is parameterically of order l/a s (a _1 ). Since a -1 ~ 1 GeV, we have not included this factor 
in our estimate. 

Our equations depend on the T mass, and make no reference to the wave function 
of the T. In fact, there are much larger nonperturbative contributions to the T(1S) wave 
function, which nevertheless, do not affect the energy of the state pTI . The reason is that the 
energy can be obtained in terms of the matrix element of the NRQCD effective Hamiltonian, 
whereas the wave function has no simple operator expansion. The situation is similar to 
other well-studied processes, such as the cross-section for e + e~ — > hadrons, where the total 
cross-section has an OPE description and much smaller nonperturbative corrections than 
other non-OPE observables such as event shapes ||33|| . 

From the results presented in Ref. [|13[], one can see three indications that the nonpertur- 
bative corrections to Eq. (f|) may be smaller than typical expectations based on Eq. (pT|). 
First, the sign of (a s G 2 ) is positive, and a large value for it increases \V c b\ as extracted 
from the inclusive semileptonic B width. Thus, for a very large nonperturbative correction 
to rrir (say, above 500 MeV or so), the consistency of the inclusive and exclusive \V c b\ de- 
terminations is significantly worse. Second, in Ref. [IB| we determined \V c b\ using both the 
mass of the Y(IS') and T(2S'). Eq. ([H]) shows that nonperturbative corrections to the n th 
5*- wave state scale as n 6 . If we assume that Eq. (^Tj) holds at least approximately also for 
the size of nonperturbative effects for higher radial excitations of the T, the nonperturbative 
contributions would be 64 times larger for the 25 than for the IS*. The fact that the resulting 
values of |V^| coincide at the 5% level at order e 2 |13j indicates that there cannot be a large 
nonperturbative correction to the mass of the T(IS'). Third, the "most nonperturbative" 
that an T can be is for it to look like a B — B bound state, which occurs just at threshold 
for decay into two mesons. The nonperturbative contribution to the B + B mass is 2A, and 
is of order 1 GeV. Assuming again that this is approximately the nonperturbative contribu- 
tion for n = 3, one finds that the Y(IS') state has a nonperturbative contribution of order 
1 GeV/3 6 ~ 1.4 MeV, an absurdly small value. 
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Finally, note that one cannot separate perturbative and nonperturbative contributions 
to a physical quantity, since the perturbation series is a divergent series.^] One can treat the 
perturbation series as an asymptotic series, and sum the terms until they no longer decrease 
in magnitude. The resulting difference from the complete answer (which is not known) can 
be called the nonperturbative contribution. If the perturbation series is poorly behaved, so 
that one can keep only a few terms in the asymptotic expansion, then one would expect 
that the error in the asymptotic expansion is large — i.e., that nonperturbative effects are 
also large. To make this more precise, consider the e + e~ total cross-section into hadrons 
at s = Q 2 , in the limit of massless quarks. By dimensional analysis, the cross-section must 
have the form 

Q 2 *(s = Q 2 ) = f MQ)) , (22) 

The dimensionless function f(ot s (Q)) includes both perturbative and non-perturbative ef- 
fects, and has an asymptotic series expansion in a s . For an asymptotic expansion, the best 
approximation to / is given by summing the series till the terms start to increase. The error 
at the minimum term, which is typically of the form 

e -c/a s (Q) „ (23) 

can be considered the nonperturbative contribution to /. 

The terms in the a s expansion for the T mass in Eq. (^) start increasing already at order 
e 2 . Using the order e term as an estimate of the error in the asymptotic series gives 100 MeV, 
which can be taken as an estimate for the size of nonperturbative contributions. When one 
combines mj with B decay rates, the perturbation series is much better behaved, and one 
expects a much smaller nonperturbative contribution to B decay rates expressed in terms of 
mx- In other words, the dominant part of the large nonperturbative contributions to mx are 
closely related to the Aqcd uncertainty in the pole mass. What is relevant for our analysis 
is the nonperturbative contributions to mx after these pole mass effects have been removed. 
The residual contributions are expected to be much smaller, of order ~ AQ CD /(m b a s ) 3 , and 
formally of the same order as 1/mf corrections which have been neglected in the inclusive 
decay rates. It is these residual nonperturbative effects for which we have used the estimate 
of 100 MeV. 

To summarize, we think that the size of nonperturbative corrections to mx relevant for 
our analysis are hardly known at the order of magnitude level (the conventional estimates 
yield hundreds of MeV, but our arguments above indicate that they could equally be tens of 
MeV only). The best way of determining them is directly from the experimental data. We 
would like to find B decay observables sensitive to mx for which parameter-free predictions 



This separation is not possible even if an infrared cutoff is introduced, since the perturbation 
series is still divergent. The n\ growth in the coefficients due to infared renormalons is removed, 
but there are other sources of factorial growth in the coefficients, such as the number of diagrams 
at n th order [p^ ]. 
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can be made. The best candidates are shape variables in inclusive decays or form factor 
ratios in exclusive decays. The reason is that observables like total inclusive rates are only 
moderately sensitive to my, and a possible discrepancy between experiment and theory 
could always be interpreted as a different value for a CKM angle, as opposed to constraining 
nonperturbative corrections to mr- The most promising observables are the ones which are 
the most sensitive to A: the mean photon energy in B —>■ Xg'j, and the form factor ratio 
Ri(w) = hy{w) /h Al {w) in exclusive B — > D*eu decay. These are discussed in Sees. V.B and 
IV, respectively. 



IV. EXCLUSIVE DECAYS 

Exclusive semileptonic B — > D*ev decays depend on four form factors, 



{D*(v',e)\ A» \B(v)) 
y/m D * m B 



h Al (w + 1) e*^ - {h M v* + h As v'») (e* • v) , (24) 



where w = v ■ v ' . In the infinite mass limit h A2 = and the other three are equal to the 
Isgur-Wise function |7j . One linear combination of the form factors is not measurable when 



the lepton masses are neglected. It is conventional J33| to define two ratios of the three 
measurable form factors 

h v {w) h Az {w) + (m D ./m B )h A2 (w) 

Ri{w) = 7 , r , Ra{w) = . 25 

h Al (w) h Al (w) 

Here we shall concentrate on the prediction for -Ri(l) at order (6 2 )blm in the upsilon expan- 
sion. Using the notation and results of Ref. [j35H , the theoretical prediction for R\ is 



where the elipses denote the known order a s /m c ,b and unknown order l/m 2 cb corrections. If 



the QCD sum rule prediction 77(1) = 0.6 ± 0.2 J3(J is approximately correct, then the A/m^ 
correction is much smaller than the A/m c term. One can avoid relying on a theoretical 
estimate for ^(1), since the deviation of R2 from unity measures i](w) directly if the \2 
subleading Isgur-Wise function is neglected. Since xi parameterizes l/m c ^ corrections due 
to insertion of the chromomagnetic operator in the Lagrangian, it is expected to be small 



and a QCD sum rule calculation supports this [0. Neglecting X2 seems to be less model 



dependent than relying on the previously mentioned result for 77(1). The relation between 
i?2(l) and r/(l) when \2 is neglected is 



"W^-^Aifrw (27) 



where z = m c /rrib. 
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The perturbative correction Cf = tja is known to order a 2 /?o |38 



CZ = VA =l + ¥- C F [4>(z) - 2] + ( ^V/9bC7 



4tt 



where a s = a s (^/m c m b ) ~ 0.27 and 



An J 



-3^-^ln(z) -6. 
1 — 2 



From results given in Ref. [|20|, we can obtain the coefficient Ci to order a 2 /3o 



C7 1 = l + gc7^) + 2] + (g)Vi 



13 j., \ 25 



(28) 



(29) 



(30) 



To combine these results with Eq. ([|), it seems simplest to reexpress a s in terms of a a (m&). 
This yields 



Thus we obtain at order (e 2 )eLM in the upsilon expansion 

fli(l) = 1.21 + 0.07e + 0.01 BL Me 2 + O.OlAx/GeV 2 + [0.31^(1) - 0.25] 



(31) 



(32) 



The quantity Ai appears in this result because we used Eq. @ to eliminate the dependence 
on the charm quark mass. The last term in square brackets is the A/m& contribution in 
Eq. (p6|), and we have assumed that it is small, so that terms of order e or Ai from this 
contribution have been omitted. A 100 MeV nonperturbative contribution to mx changes 
the prediction for Ri by ±0.021, which is about twice the order (e 2 )eLM correction. Probably 



the largest uncertainty is due to the neglected l/m cb corrections p9[| , which need to be better 
understood before Eq. (|3"2"|) may be trusted at the percent level. 

One can compare Eq. (|32|) with the conventional treatment. The size of the order a 2 s f3o 
terms is usually quoted in terms of the BLM scale. The BLM scale of C\ = t)a is Q.bl^/mjnl 
and that of C\ + C 2 + C 3 = i] V is 0.92^m c m b [PS| . Surprisingly, the BLM scale of C\ is much 
smaller, 0.16yjm c m b , and that of G\jG\ in Eq. Q3lD is only slightly larger, 0.13m b . These 
low BLM scales show that the conventional perturbation series for R\ is poorly behaved. 
The perturbation series in Eq. (|32|) cannot be converted into a BLM scale since it mixes 
different orders in a s . However, to compare with conventional results one can define an 
effective BLM scale by using the ratio of the (c 2 )blm and e terms, just as the conventional 
treatment used the ratio of the (o; 2 )blm and a s terms. The effective BLM scale of Eq. ([32]) 
is 0.89m b , which is large and allows for a reliable expansion. 

The present experimental results are [HO] 



Rt = 1.18 ±0.30 ±0.12, 



Ri 



0.71 ±0.22 ±0.07. 



(33) 



This is consistent with Eq. (03), although the errors are sizable. If R^^iX) wm be measured 
in the future with one (or at most a few) percent errors, then it may be possible to constrain 
nonperturbative contributions to m T at the 100 MeV level, which would be very interesting. 
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To summarize, Eq. (|32|) predicts the form factor ratio R\ independent of A and the quark 
masses. If this is tested and works at the few percent level, then we can make better and 
still model independent predictions for B decay form factors in general, including the ones 
to the orbitally excited mesons Di and D\ |4l]| . Even if nonperturbative contributions to 



mx cannot be constrained too well, our results show that expressing B decay rates in terms 
of mx is a useful and theoretically consistent parameterization. 

V. INCLUSIVE B DECAY RATES 

In this section, we apply the upsilon expansion to inclusive B decays. The h quark de- 
cay rate has long been known for arbitrary masses of the four fermions participating in the 
weak decay, both at tree level and at order a s The order Aq CD / m 2 , nonperturba- 
tive corrections will be taken into account for all decay modes, however, a s corrections to 
these will be neglected. The numerical results below correspond to the choice fi = irtb and 
a s (mh) = 0.22, and //-dependence will be discussed as well. 

Nonleptonic decays are more complicated than semileptonic decays. It is less clear 



whether local duality holds in nonleptonic decays, since there is no smearing variable |44| , f45 
There are also corrections not present in semileptonic decay because the four-quark operators 
run in the effective theory below the weak scale. While the complete series of leading |46| 



and subleading [47[ logarithms have been summed, we cannot consistently use these results, 
since that would involve summing series of the form [a s \ia(mw / 'mb)} n and a s [a s \n{mw /mb)] n 
(where ln(mw/mb) ~ 2.8) without summing q.^ +1 [3q (where (3q ~ 9). Instead, we follow an 
approach somewhat similar to Ref. based on the observation that the largest correc- 
tions arise at low orders in a s . This is quite reasonable, as the order [a s \n{mw / rnb)] 2 term 
dominates the sum of leading logs (there is no order a s \n(mw jwib) term due to color con- 
servation), and the order a s term dominates the sum of next-to-leading logs. For b — > c 
nonleptonic decays we include the order a s correction, and at order a 2 the terms containing 
[\ia(mw / mb)} 2 , \n(mw / iTib) , and (3 .f\ For b — > u nonleptonic decays we include the order a s 
and order [a s \n(my/ / nib)} 2 corrections only. 

A. Semileptonic B decays 



Implications of our method for semileptonic B decays were studied in Ref. [|TI| , and here 
we repeat the results briefly for completeness. Estimates of the small non-BLM contributions 
to the decay rates at order a 2 , which were taken into account in Ref. [13||, will be neglected 



here. For massless leptons, the a s correction to free quark decay is known analytically [|50 
and the order a 2 /3o correction has been determined in Ref. Jl 1 



^It was argued that other order a 2 terms might also be unexpectedly large due to Coulomb 



enhancement [49|. 
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1. B -> X u eu 



The B — ► X u ez/ decay rate at order (e 2 )eLM in the upsilon expansion is [ff 



1927T 3 



0.115e- 0.035 BL Me 
9A 2 - Ax 



2(m T /2)< 



+ 



(34) 



The perturbation series, 1 — 0.1 15e — 0.035blmc 2 , is far better behaved than the series in 
terms of the pole mass, 1 — 0.17e — 0.13blmc 2 , or the series expressed in terms of the MS 
mass at /i = m&, 1 + 0.30e + 0.19BLMC 2 - The uncertainty in the B decay rate using Eq. ( |34|) 
is much smaller than that in Eq. (|I|), both because the perturbation series is better behaved, 
and because the T mass is better known (and better defined) than the b quark mass. 

Eq. (|34D yields a relation between \ V u b\ and the total semileptonic B — > X u ev decay rate 
with very small uncertainty,^ 



\V uh \ = (3.06 ± 0.08 ± 0.08) x 10" 
(B{B -> X u ev) 1.6 ps' 



1/2 



0.001 



tb 



(35) 



where we have used A 2 = 0.12GeV 2 and Ai = (—0.25 ± 0.25) GeV 2 . The first error is 
obtained by assigning an uncertainty in Eq. (34) equal to the value of the e 2 term and the 
second is from assuming a 100 MeV uncertainty in Eq. (|]). The scale dependence of |V^&| 
due to varying fi in the range m^/2 < \i < 2mb is less than 1%. The uncertainty in Ai 
makes a negligible contribution to the total error. It is unlikely that B(B — > X u ev) will 
be measured without significant experimental cuts, for example, on the hadronic invariant 
mass [51|. Our method should reduce the uncertainties in such analyses as well. 



2. B -> X c e v 

The B — ► X c ev decay rate at order (6 2 )blm is [J3 
T(B 



X c ev) 



G 2 F \V cb \ 2 fm r \ 5 
192tt 3 V 2 J 



0.533 



1 - 0.096e - 0.029 BL Me 2 
- (0.28A 2 + 0.12Ai)/GeV 2 



(36) 



For comparison, the perturbation series in this relation when written in terms of the pole 
mass is 1 — 0.12e — 0.07blmc 2 - Equation (j36l) implies 



\V cb \ = (41.6 ± 0.8 ± 0.7 ± 0.5) x 10" 3 

f B{B^X c ev) 1.6ps\ 1/2 



x 



^QED 



0.105 



(37) 



9 In Eqs. (35) and (|37D, we have included the non-BLM terms of order e 2 as well (see Ref. |T 
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where t/qed ~ 1.007 is the electromagnetic radiative correction. The uncertainties come 
from assuming an error in Eq. (|36|) equal to the e 2 term, the 0.25 GeV 2 error in Ai, and 
a 100 MeV error in Eq. (f|), respectively. The second uncertainty is reduced to ±0.3 by 
extracting Ai from the electron spectrum in B — > X c ei>] see Eq. (p8|). The agreement of 
|V c fe| with other determinations (such as exclusive decays) is a check that nonperturbative 
corrections to Eq. (|4]) are indeed small. 
In Ref. 



A and Ai were extracted from the lepton spectrum in B — >• X c ev decay. 
With our approach, there is no dependence on A, so we can determine Ai directly with small 
uncertainty. Considering the observable R\ = J 15GeV E e (dT/dE e )dE e / J 15GeY (dT/dE e )dE e , 
a fit to the same data yields 



Ax = (-0.27 ± 0.10 ± 0.04) GeV 2 . 



(38) 



The central value includes corrections of order a 2 /?o [p3" |. The first error is dominated by 
1/m 3 corrections |j4]. We varied the dimension-six matrix elements between ±(0.5 GeV) 3 , 
and combined their coefficients in quadrature in the error estimate. The second error is from 
assuming a 100 MeV uncertainty in Eq. (£|). The central value of Ai at tree level or at order 
a s is within 0.03 GeV 2 of the one in Eq. Q3"8|). This value of Ai can be used consistently in 
theoretical expressions up to order a 2 /5 |5 



3. B — > X r tv,B^ X n t v and R T 



The order AQ CD /m 2 corrections are taken from [ 56] , |57 1 , and the order a 2 /?o corrections 
from |j48"l . The B — > X c rv decay rate at order (€ 2 )blm is 



T(B 



X C TV) 



192tt 3 V 2 J 



0.070e - 0.016 BL Me 2 
(0.59A 2 ±0.25Ai)/GeV 2 



(39) 



For comparison, the perturbation series in this relation when written in terms of the pole 
mass is 1 - 0.097e - 0.064 BL Me 2 . 

For the B — > X u rv decay rate we find 



T(B -> X u tu) 



G 2 F \V ub \ 2 ( m r \ l 
192tt 3 I 2 J 



0.361 



1 - 0.08e - (0.34A 2 - 0.02A!)/GeV s 



(40) 



For comparison, the order e correction written in terms of the pole mass is 1 — 0.16e. 

The ratio R T = B{B — > X c ru) /B(B — > X c eu) is interesting because it directly probes 
lepton mass effects, and it can be predicted precisely in the standard model. The largest 
uncertainty in the prediction is Ai, and varying it in the range given in Eq. (|38f) yields 0.220 < 
R T < 0.227.Q Thus, in the standard model, measurements of R T can be used to constrain 
Ai H,||,||. The most recent ALEPH data, B(B -> X c tu) = (2.72 ± 0.20 ± 0.27)% ||, 



10 



For this calculation we have set B(B 



X c ev) 



10.5%, and neglected its uncertainty. 
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implies at the la level Ai < — 0.36GeV 2 . However, at the 1.5a level this upper bound 
already becomes positive, which is not very interesting since Ai is expected to be negative. 
R T is also interesting for constraining new physics, e.g., possible (tree level) charged Higgs 
contributions in the large tan (3 region pT . 



B. B - X sl 

The mean photon energy over the region Eq < E 1 < _E™ ax in inclusive B — > X s 7 decay 
for sufficiently small Eq is sensitive to A and independent of Ai^ when the decay rate is 
expressed in terms of and A instead of m;, f62fl . We expect that stringent constraints can 



be obtained on the nonperturbative corrections to mx when a measurement of the photon 
energy spectrum is available with E lowered to around 2 GeV p2[|63 |. 



Our method reduces the theoretical uncertainties in the total B — > Xgj decay rate only 
moderately. The B — > X s ^ decay rate is usually normalized to the semileptonic B — > X c ev 
rate to eliminate the m\ factor from the prediction. The two biggest ingredients of the total 
theoretical error are the scale dependence of the B — > X c ev rate computation itself (without 
the m\ factor), and to a lesser extent the normalization to the B — > X c eu rate due to the 
uncertainty of m c /mb and the experimental error. 

C. Nonleptonic B decays 

The order AQ CD /m 2 corrections are taken from ||||. We neglect the strange quark mass, 
and the deviation of |V^ S | 2 + \V c d\ 2 and \V u d\ 2 + \ V US \ 2 from unity. 
The most interesting decay mode is 6 — > cc(s + d). 

T(B -> X c5[s+d) ) = °l^f (^) 5 0.212[l + 0.156e + (0.074 B lm + 0.157L 2 + 0.046L)e 2 

- (0.45A 2 + 0.41Ai)/GeV 2 ] , (41) 

The variable L — 1 is introduced to distinguish between terms that come from different 
powers of ln(mvp/w&). Equation ( |4"T{ ) has been written with the renormalization group 
scaling factor expanded in powers of a s ln(M^/m&). One can effectively include the complete 
leading logarithmic series by replacing the 0.157L 2 term by 0.113L 2 . This replacement 
also holds for Eqs. (f42|), fl43|) and (|44|) . The upsilon expansion does not affect the order 
e 2 L and e 2 L 2 terms, i.e., they are the same as in the conventional a s expansion. The 
behavior of the order e and (e 2 )eLM terms have improved significantly. When written in 
terms of the pole mass, these terms are 1 + 0.199e + 0.151blm£ 2 {HI- The cancellation we 



find is rather nontrivial. The expansion of ml when reexpressed in terms of mx generates 
positive contributions at order e 1,2 . This is why there was a significant cancellation in the 
perturbation series for b —>■ uev decay, since the conventional perturbation series expansion 
for this decay had negative e 1,2 terms. On the other hand, the conventional perturbation 
series for nonleptonic b — > cc(s + d) has positive e 1,2 terms, and reexpressing the overall 
ml factor in terms of mj- only makes them larger. It is the expansions of the tree level 
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b — > cc(s + d) phase space, which generates large negative contributions at order e 1 ' 2 to make 
the expansion reasonably well-behaved. 

For nonleptonic b — > cu(s + d) decay, the rate is 



T(B - X, 



cu(s+a 



G 2 F \V cb \ 2 fm r \' 
64tt 3 V 2 J 



0.533 



1 - 0.026e + (-0.006 B lm + 0.157L 2 + 0.099L)e 2 



- (0.28A 2 + 0.12A!)/GeV 5 



(42) 



When written in terms of the pole mass, the order e and (e 2 )sLM terms are 1 — 0.048e — 

0.045BLM6 2 . 

The decays mediated by b — > u transitions are much less important, and therefore we 
keep fewer terms than before. 



r(B - x v 



G 2 F \V ub \ 2 fm r V 



1 - 0.045e + 0.157Z/e 



2 2 



^uu{s+d)j 647r 3 y 2 J 
For comparison, the order e correction written in terms of the pole mass is 1 — 0.099e. 



uc{s+d) j 



G 2 F \V uh \ 2 fm T V 
64tt 3 V 2 J 



0.533 



1 + 0.116 + 0.157^ 



2,2 



(43) 



(44) 



This is the only case where the upsilon expansion increases the order e correction compared 
to when it is written in terms of the pole mass, 1 + 0.09e. 



D. Semileptonic branching ratio and charm counting 

The theoretical prediction for the semileptonic B branching ratio Bsl was discussed 



in ||64|| . It has been emphasized repeatedly that the calculation of the b —> cc(s + d) rate is 
the least reliable. With our expansion in Eq. (|4~ID , the ratio of the (e 2 )eLM correction to the 
order e term is significantly smaller than that calculated in terms of the pole masses, so the 
b — > cc(s + d) rate should be under better control. With our central values for the inclusive 
decay rates given in the previous section we obtain 

23 SL = 11.8%, n c = 1.22, (45) 

for the semileptonic branching ratio and the average charm multiplicty, which is consistent 
with other recent analyses. Here we have used \V u b/V c b\ = 0.1. The sensitivity of these 
results to Ai, for example, is small when it is varied in the range obtained in Eq. d38|) . If 
we had not expanded the series of leading logarithms, n c would remain unchanged, but Bsl 
would increase to 12.1%. Including a +100 MeV nonperturbative contribution to the T 
mass changes these values by a negligible amount. 



VI. CONCLUSIONS 

We have shown that inclusive and exclusive B decay rates can be predicted in terms of 
the T(15) mass instead of the b quark mass. It is crucial to our analysis to use the upsilon 
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expansion in e rather than the conventional expansion in powers of a s . Our formulae relate 
only physical quantities to one another. They result in smaller theoretical uncertainties 
than existing numerical predictions, and the behavior of the perturbation series is improved. 
Moreover, the uncertainties can be estimated without resorting to cumbersome arguments, 
and they can be checked using the experimental data. 



In a previous paper |13| , we had shown how our method works for inclusive semileptonic 
B decays. In this article, we have shown how it applies to exclusive decays and to nonleptonic 
decays. In the new cases studied, the upsilon expansion also improves the behavior of the 
perturbation series. The biggest theoretical uncertainty at present is a reliable bound on 
the effect of nonperturbative contributions to the T mass (other than the linearly infrared 
sensitive piece of the pole mass). The best way of determining this seems to be directly 
from the experimental data. The comparison of V c b extracted using the upsilon expansion 
with previous determinations already suggests that these contributions are not large. We 
are investigating whether other processes can be used to provide an independent test of this 
result. 
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